New mechanism for non-trivial intra-molecular vibrational dynamics. 
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We investigate the time evolution process of one selected (initially prepared by optical pumping) 
vibrational molecular state S, coupled to all other intra-molecular vibrational states R of the same 
molecule, and also to its environment Q. Molecular states forming the first reservoir R are charac- 
terised by a discrete dense spectrum, whereas the environment reservoir Q states form a continuous 
spectrum. Assuming the equidistant reservoir R states we find the exact analytical solution of the 
quantum dynamic equations. S - Q and R - Q couplings yield to spontaneous decay of the S and R 
states, whereas S - R exchange leads to recurrence cycles and Loschmidt echo at frequencies of S - 
R transitions and double resonances at the interlevel reservoir R transitions. Due to these couplings 
the system S time evolution is not reduced to a simple exponential relaxation. We predict various 
regimes of the system S dynamics, ranging from exponential decay to irregular damped oscillations. 
Namely, we show that there are four possible dynamic regimes of the evolution: (i) - independent of 
the environment Q exponential decay suppressing backward R -S transitions, (ii) Loschmidt echo 
I , regime, (iii) - incoherent dynamics with multicomponent Loschmidt echo, when the system state ex- 

changes its energy with many states of the reservoir, (iv) - cycle mixing regime, when the long term 
system dynamics appear to be random. We suggest applications of our results for interpretation of 
femtosecond vibration spectra of large molecules and nano-systems. 
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q ; I. INTRODUCTION 

A quasi-stationary state in quantum mechanics occurs as a result of overlapping (superposition) of an initially 
qq ' localised (stationary) state with the states (i.e., wave functions) from continuous spectrum formed by either free 
£N) \ states of the same system (separated from the initial state by a potential barrier, as in the case of a-decay), or by the 

• system environment (see classical papers [2, 0). Following this, common wisdom ascribes irreversible evolution of 
quasi-stationary states by coupling the states to a reservoir with continuous spectrum 1 , 1 ■ With a model reservoir 

. formed by a sea of harmonic oscillators, this approach is at the heart of the theory of quantum dissipative systems 

o: 1, 1, 0- 
on- 

An opposite limit is considered in the theory of transition states which is widely used to treat various chemical 
dynamics problems 1,1- In this microscopic approach one has to choose properly a set of internal degrees of 
• i-H , freedom forming a so-called reaction path, and a small number of transverse degrees of freedom which are coupled 

• to the longitudinal reaction coordinate. This microscopic description is feasible in practice up to a few dozens of the 
transversal degrees of freedom. For larger systems the microscopic approach becomes useless (since even with modern 
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computer powers it is almost impossible to solve a system of dynamic equations with 10 4 transversal degrees of freedom 
for all eigen states to restore multidimensional potential energy surfaces). However many nano-systems with 10 2 — 10 4 
degrees of freedom, interesting from their practical importance and the associated theoretical challenges, belong to an 
intermediate case when both of the above mentioned approaches (macroscopic theory of quantum dissipative systems 
and microscopic theory of chemical dynamics) do not work. Evidently to cover these very complex phenomena, it 
is necessary as a first step at least to choose an appropriate simple model. One approach is to borrow concepts 
from other physical systems with dense discrete spectra, for example nuclei. It was shown in [l(J, [ll[, [l3|, 
[Til ], that statistical description of such systems does not require any detailed information about its spectrum, but 
only a few universal spectral characteristics (like interlevel spacing distribution function), determined by a random 
Hamiltonian matrix. This approach is a very convenient tool to describe spectral chaos and many other global features 
of the behavior, but it says almost nothing about quantum dynamics, in which we are interested in this paper. Our 
motivation is not a pure curiosity. As a matter of fact quantum dynamics of various systems (ranging from relatively 
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small molecules in apre-dissociation condition [Tj|, [T|| up to lar ge p hotochromic molecules and protein complexes 
[13 > [3) [13]; IHI) [HI)) or molecules confined near interfaces |23j] (see also (22j) is an active area of experimental 
research. Femtosecond spectroscopy data (which allows the study of the time evolution of one initially prepared by 
the optical pumping state) manifest a variety of possible dynamic regimes including not only weakly damped more or 
less regular oscillations but also very irregular long term behavior with a number of peaks corresponding to a partial 
recovery of the initial state population. 

Seemingly irregular damped oscillation regimes observed in such systems cannot be explained theoretically in the 
frameworks of widely used models with reservoirs possessing continuous spectra [24[ , [j| , Q , Q , [H • Indeed in the 
case of a system coupled to the continuous spectrum reservoir, only a smooth crossover between coherent oscillations 
and an exponential decay is possible upon increasing of the coupling. Nevertheless, as it was shown in the papers 
cited above, generic complex dynamics is observed in the systems with characteristic inter-level spacing of the order 
of 10 cm -1 , when the measurements are performed in the range of sub-picoseconds, or femtoseconds. To explain these 
weakly damped oscillations semi-empirical models have been proposed p2l |. pEl ]. assuming more or less arbitrarily 
that the system interacts not only with the environment (possessing continuous spectrum in the agreement with a 
common belief), but also with a few weakly damped discrete vibrational levels. However, these models providing a 
possible mechanism for weakly damped oscillations, do not explain irregular, random-like dynamic evolution. It is 
worth noting again that a generic feature of systems with such irregular behavior is the existence of the dense but 
discrete vibrational spectra, about 10 2 — 10 3 levels, and characteristic inter-level spacing of the order of 10 cmT 1 . 
This generic feature is present in all systems with complex and irregular vibrational relaxation. Motivated by these 
observations, our intent here is to propose a simple (but yet non-trivial) model of a system coupled to a reservoir 
with a discrete spectrum, and to examine joint system-reservoir evolution, i.e., recurrence cycles, when the energy is 
flowing back from the reservoir to the system, and to apply this model to intra-molecular vibrational relaxation. 

Vibrational relaxation in large molecules is one of the most relevant processes in chemical dynamics. The system 
time evolution includes as its step intra-molecular energy transfer via transitions from a specially prepared and selected 
initial molecular state (system S) to a large finite number of distinct levels of the same molecules. This intra-molecular 
evolution should be supplemented, of course, by the interaction of all these states that posseess a continuous spectrum 
environment, which ultimately yields to conventional exponential relaxation [26[ - [29(. To describe theoretically these 
radiationless transitions, usually one has to rely on approximate treatment of the system - reservoir interactions, e.g., 
on random phase approximation. In this case the system time evolution is determined by the decay rate constant, 
which is calculated accordingly with the famous Fermi Golden Rule. Unfortunately such an approach is valid only in 
the limit of a relatively weak system - reservoir coupling, when calculating perturbatively transition probability, one 
may neglect shift of the energy levels [3(| ■ 

In our paper we develop another approach to analyse theoretically the problem of intra-molecular vibrational 
dynamics. Instead of a perturbative solution to the quantum dynamic equations for a system coupled to an arbitrary 
reservoir, we find the exact solution for a simple model of a reservoir spectrum (i.e., a set of final states). The 
model keeps the main essential feature of intra-molecular vibrational dynamics, namely that realistic reservoir should 
possess discrete and dense spectrum of its states. Surprisingly enough, scanning the literature we could not find any 
paper treating theoretically such a model. We do believe that the basic ideas inspiring our work can be applied to a 
large variety of interesting nano-systems. In the recent short publication [3l| such a model of a system coupled to a 
reservoir with dense discrete spectrum was proposed, and under assumptions put forward by Zwanzig [32j (equidistant 
spectrum of the reservoir and system-reservoir coupling independent of reservoir state quantum numbers) its exact 
analytic solution was found. Here we generalise this approach to rationalise intra-molecular vibrational relaxation in 
terms of time evolution of quasi-stationary states. The key feature, indispensable to provide energy exchange between 
the system and reservoir, is backward energy flow, which is ignored in all models with continuous spectra of final 
states @ . These backward transitions (and corresponding recurrence cycles or Loschmidt echo phenomenon) are the 
main new physical ingredients of our approach, and they are responsible for non-trivial time evolution. For a fairly 
broad class of molecules, such a non-trivial time evolution indeed was observed experimentally (see the monographs 



in [33[ , [34[ and references therein) by modern femtosecond spectroscopy technique. Our model seems to be adequate 
to rationalise these experimental observations, in a sense that we are not aiming at the best quantitative agreement 
of theoretical results and experimental data but at the identification, understanding and explanation of important 
sources of non-exponential time evolution. 
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II. BASIC THEORETICAL NOTIONS FOR A MODEL WITH TWO RESERVOIRS. 



Our analysis is based on a simple but very general observation known from standard quantum mechanics. For a 
selected unperturbed energy level e° (in what follows we will term the level as a system) coupled to a discrete spectrum 
reservoir (with its unperturbed energy states {e°}), the Hamiltonian matrix contains besides the main diagonal, only 
one row and one line of non-zero matrix elements: 
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Therefore the corresponding secular equation (its roots determine the coupled system-reservoir eigenvalues) has the 
following deceptively simple form [3l|, (3f| 

^) = ^-EtS^ ' (2) 

n n 

where we count the energy levels from e^, and to get such a compact form for the function F(e) we use the orthogonal 
basis of the reservoir states, i.e., all matrix elements between the reservoir states are zero. With the same approach 
we can include all ingredients of radiationless molecular transitions into a simple model, namely 

• system S which is selected molecular level 

• discrete dense spectrum reservoir R, which includes all other than S molecular states 

• S — R coupling characterizing by the coupling matrix elements {C„} 

• second reservoir (environment) Q with continuous spectrum, which leads to a decay of the system S (with its 
rate constant T s ) and of the reservoir R states (with corresponding constants {T n }). 



One can derive formally an exact secular equation similar to the ((2]) for the system S and reservoir R eigenvalues, which 
additionally includes the environmental level broadening (due to S — Q, and R — Q couplings). The corresponding 
determinant contains (besides the diagonal matrix elements E + iT Sl and {E — e„ + iT n }) only one non-zero row and 
one line with non-zero S — R coupling matrix elements. The corresponding secular equation reads as 

F(E)^E + iT s -^ E _^ +iTn =0. (3) 
Of course for arbitrary functions C„, £° , and r„ , the sum in the Eq. ([3]) cannot be calculated analytically. To proceed 



further on we are following the Zwanzig idea [321 



(already broadly used in the theory of radiationless transitions [2f 



29j ] . and time resolved spectroscopy [34]] > [3q|. 37l]) assuming that 

e ° =n;C n = C;T n =T, (4) 

where as above, we count all energy levels from the system S energy unperturbed by its coupling to the reservoir 
levels, and the unperturbed reservoir R interlevel spacing is chosen as the energy unit. It is worth noting also that 
the approximation of equidistant R reservoir spectrum is not a completely artificial one. Indeed, for any system with 
a sufficiently large finite number of degrees of freedom, mutual level repulsion unavoidably favours to more or less 
equidistant interlevel spacing. A more realistic model will not affect our qualitative conclusions, and transparency is 
worth a few simplifications. If the Zwanzig model assumptions are granted, the secular equation ([3]) can be written 
in the following compact and easily solvable form 



F(E) =E + iT s - ttC 2 cot[7r(£ + iT)} = . 



(5) 
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System S and reservoir R energy states can be regarded as quasi-stationary ones if [38| 

T s < 1 ; T < 1 , (6) 
then complex solutions to the secular equation ([5|) are 

E n =e n - Z7„ , (7) 



where 



and 



^2 \ I /t~i f n2 en cot(7re„) + ttC ] cot(7re») 

e„ = ttG cot(7re„j + (1 s - 1 j -; r— — , (8) 

e„ cot(7re„) + 7rC/ 2 + (l/7r) 2 



7n = r+(r fl -T) 57 ~i _i_ n2 ~j~7T~/ \ ' ( 9 ) 

£„ COt(7T£ n ) + 7rC 2 + (1/7T) 

Note that the shift of the zero level energy corresponds to a change of variables in a rotating frame, and it does 
not affect system dynamics in which we are interested in this work. By a simple inspection of the Eqs. 1(5)) - © we 
conclude that S — Q and R — Q couplings displace the system S and reservoir R levels into the lower complex semiplane 
with energy dependent decay rates. S — Q, R — Q, and S — R interactions yield to variations of all energy levels. It 
is easy to see that the energy level shifts "e„ — n" are mainly due to S — R coupling. It means that independent of 
the coupling strength C there is one energy level in each interval [n , n + 1] . 



III. SYSTEM AMPLITUDE 



To analyse the system time evolution we have to solve the time dependent equations of motion (i.e., the corre- 
sponding Heisenberg equations of motion) for the system state amplitude a s . Our Hamiltonian model describing 
the system S and the discrete spectrum reservoir R (while producing decay of the quasi-stationary states couplings 
with the continuous spectrum reservoir Q are taking into account by introducing the decay rates T s and r n into the 
equations) 

H = 6°6+6 s + ^ b nbn + Cn(btbn + b s b+) , (10) 

n n 

where bj, are corresponding creation operators. Time dependent wave functions ^ s (t) of the Hamiltonian (JTUJ) 
can be expanded over the unperturbed (uncoupled) eigenfunctions of the system $ s and of the reservoir states 

*,(t) =a s (t)<t> s +Y,a n (t)<S> n (11) 

n 

with time dependent amplitudes a s (t) , a n (t) . These time dependent amplitudes satisfy the corresponding Heisenberg 
equations of motion. Combining everything together we end up with a general case (i.e., not applying immediately 
the Zwanzig approximation), the equations read (in units with % = 1, and dots steam for time derivatives) [3l| . [35| 

ids = -«r s a s + ^2 C n a n ; ia n = (e° - ir n )a n -I- C n a s . (12) 

n 

The solution to the Eq. p2[) supplemented by the initial conditions 

a s (0) = l;a„(0) = 0, (13) 

can be formally found as 

i ^-oo / pi \ 1 
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The amplitude can also be expressed as a sum of residues in the poles of the integrand in the Eq. (|14p which are the 
roots of the secular equation 



^ = y-, exp(-iEt) 



dF/dE 



(15) 

F(E)=0 



Let us apply now the Zwanzig approximation ((4]). In this case the formal solution for the system amplitude (|15[) can 
be written down in the explicit analytical form 

a s(t) = i + 7r 2C2 X + C -2f l R +iT \2 = 2 E e M-lnt)[B n cos(e„i) + S n sin(e n t)] , (16) 

n— — oo ^ s ' n— 

where 

£?„ = - (r, - 7 „) 2 + C 2 (l + irC 2 )} ; , (17) 

and 

S n = 2Be n {T s - ln ). (18) 

The coefficient £? in these formulas reads as 

B = C\el (r s - 7 „) 2 + C 2 (l + tt 2 C)] 2 + 4e 2 (r s - 7,,) 2 ]- 1 . (19) 

The above expressions (fT6|) - (fT9|) give the exact solution to describe quantum dynamics of the system S coupled to 
the Zwanzig reservoir R and continuous spectrum environment reservoir Q. 



IV. RECURRENCE CYCLE AMPLITUDE REPRESENTATION 



Presented in the previous section exact results are not very useful in practical terms, because of a poor convergence 
of Fourier series entering these expressions. Luckily one can transform the results into a form more convenient for 
further calculations (both analytical and numerical) by replacing the Fourier expansion with the expansion over partial 
recurrence cycle amplitudes. The method of how to do it, developed in the papers [3l|, [3!|, (3f|, is based on the 
Poisson summation formula, (see e.g., [HI). To transform a s (t) as it is given in the IT6|) . we introduce the following 
formal identity for the S function 

-l 

S(x - n) . (20) 

The Eq. ([20|) allows us to replace the discrete variables n and E n by their continuous counterparts x and E(x) with 
E(n) = E n . Thus we end up with 

o«W = (;r) exp(-iEt)S(E - E n ) = (21) 



5{E{x))=Y J 



dE 
dx 



^ J dE \dE~) \~cT) &c P\~ iE ( t ~ 2/c7r ) + i2k7T ( E - x ) 



k— — oo 

The expression f|2ip represents system evolution as a sum over recurrence cycles 



,(*)= E ^\t-2k-K), (22) 



k— — c 



where the partial cycle k amplitude (for Zwanzig reservoir R) reads as 

/OO / Tp , --p • S~l2\k — 1 

J E {E + iT a + - 2fc7r)] . (23) 
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In the initial cycle fc = 0, there are two poles E± = — i(T s ± ttC 2 ) in the integrand for ai "*, and for t > only the 
pole El contributes, yielding to the exponential relaxation law for the ai°^ (t) 

ai°)(i)=exp[-(r s +^C 2 )]. (24) 

We see that S — R coupling gives an additional (with respect to the environment Q) decay channel with its rate 

T R = nC 2 . (25) 

In a higher order recurrence cycles k > 1, the pole E*_ disappears and the residue in the (k + l)-th order pole E1_ 
leads to 

a { s k) (t - 2kn) = exp[-2fc7iT] exp[-L s (t - 2ki:)]a% ) , (26) 
where for the sake of compactness 

a#(*) = - 2rfl(t ~ 2fc7r) 4_i[2r fl (f - 2kn)}0(t - 2kv) . (27) 

Here fc > 1, L\_ l is adjoint Laguerre polynomial [4l|, and (z) is a step function defined as Q(z < 0) = 0, and 
0(z > 0) = 1. Since the pole El exists only for t > 2kir in the representation of the recurrence cycle partial 
amplitudes for a s (t) only the cycles with k ranging from to [t/27r] ([x] stands for the integer part of x) contribute 
into the sum entering equations (f2Tj) - ([22]) . 



Similarly the amplitudes for the reservoir R states can be written as 



where for fc = 



[t/2n] 

l (t)= J2 4 fc) (*-2fc7r), (28) 

fc=0 



4 0) W = I ttt (exp[-(r s +r fl )t] -exp[-(m + r)i]), (29) 



n + i(r s + t r -r) 



a^\t-2kir) = -iCexp[-(in + T)(t-2kn)] I d£ exp((m + Tjt')^^ (t') . (30) 



and for k > 1 



It is worth noting that although transformation from the Fourier representation into the partial amplitude repre- 
sentation mathematically is merely a change of variables, there is some physics behind it. By inspection of both 
representations one can see that the Fourier expansion coefficients arc determined by the real eigenvalues, while the 
partial amplitudes depend mainly on imaginary pole displacements (which in their own turn are related to the S — R 
coupling or decay rate Tr). This is not a minor technical issue but an essential physical difference, because in fact 
it means that the partial recurrence cycle amplitudes take into account effects of all R states. Physics occurs as a 
result of it. Indeed quantum interference of transitions from different R states suppresses effectively back transitions 
from the R states into the system S, and yields to a new decay channel. Its decay rate is determined by the tran- 
sition probability into all R states. Back transitions from the R states into the system lead also to the spontaneous 
Loschmidt echo at each recurrence cycle. Instead of an individual level-to-level transition picture (satisfying standard 
detailed balance relations) we have to deal with cooperative transitions from S into all R states. Tr is the rate of 
these transitions, and for the backward transitions from R to S the effective rate constant is of the order of unity, 
because the reservoir R should "wait" for synchronisation of its states to have backward R to S transition. 

This non-trivial behavior manifests itself in a fine structure of the Loschmidt echo signals. Utilising known properties 
of the Laguerre polynomials (4l| we arrive at the conclusion similar to that found in [3l| , [39l | , (35| for a single reservoir 
Zwanzig model. Namely, that in the fc-th recurrence cycle, the Loschmidt echo signal has k components, Iti with 

I = 1, , fc. The component Ikk is the most intense one. Because different R states exchange their populations 

with the system S not at the same time, the echo signal is broadened. The full width of the signal in the cycle fc is 
about 4fc. In its own turn when fc grows, the Loschmidt echo signals for the neighboring cycles start to overlap, if 
fc > fc c = 7r 2 C 2 , and for fc ^ fc c time evolution looks as irregular (quasi-stochastic) damped oscillations. 
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It is also interesting to compare the formally exact presentation of a s (t) in terms of the partial cycle amplitudes 
(f2"2")l - (j2"3"|) with its approximate treatment according to the Golden Fermi Rule. The equations of motion (fT2"]) written 
in the equivalent integral form can be rewritten as the quantum Langevin equation for kicked harmonic oscillations 
with friction 

a s (t) = -T s a s - f dt'a s {t') V C 2 n exp[-i(e£ - iT n ){t - t')] . (31) 

The terms entering this equation have transparent physical meaning: the resonance transition to the level n = of 
the reservoir leads to coherent oscillations with frequency oc C, whereas non-resonance (n ^ 0) transitions produce 
effective friction and exciting force. However, strictly speaking, there is no stochasticity in the equation, and we have 
fully deterministic energy levels and therefore time evolution. 

In the Golden Fermi Rule approximation we replace 

/oc 
p(e)de, (32) 

n -°° 

with p(e) being unperturbed (by interactions) density of states. Furthermore, if C^(e°)p(eo) is a slow varying (in 
comparison to exponential) function, the integrand in the (|31[) - Q32[) can be approximated by the mean value of this 
function. Thus, according to the Golden Fermi Rule 

a s (t) = -(T s + T GFR ), (33) 

where Tqfr = 27r(C 2 (e)p(e)), and therefore in this approximation the memory about all recurrence cycles is lost. For 
a single reservoir Zwanzig model [3^ | this approximation corresponds to keeping only one term k — in the Poisson 
summation formula. 



V. TIME EVOLUTION REGIMES 



The partial cycle amplitudes a s are represented according to the expressions (|26|) . (|27|) as a product of three 
factors. The first factor is related to the R — Q coupling. Physically it describes the Loschmidt echo signal reduction 
due to spontaneous decay of the R states in the previous (with respect to the cycle k under consideration) cycles. In 
the strong coupling limit we are interested in, i.e., for Tr 3> 1, this factor does not depend on T s . Indeed in this limit 
in the major part of the recurrence cycle only the reservoir R states are populated, whereas the system population is 
negligibly small. The second factor in the Eqs. (|26|) . (|27p is related to the system decay in the same cycle k. Naturally 
this factor does not depend on the decay rate of the reservoir R states. Finally the third factor in the (f2T)|) . (f2"T)l is 
related to the decay rate the system at T s = T = 0. 

These relations (PIBl) - ([2"T]) and their interpretation are our main results in this paper. Analysing the relations we 
end up with the following classification of the possible time evolution regimes. Namely 



• Exponential decay in the initial cycle k = 0, when the Loschmidt echo intensity is exponentially small 

r » i- . (34) 

If Tr 3> T s then the decay rate r# + T s ~ does not depend on the reservoir Q. Physically it means that in 
this limit the rate of radiationless intra- molecular transitions is independent of the molecular environment Q. 
The states of the R reservoir with \n\ < Tr gives the main contributions into the S state decay rate (dominated 
in this limit by the S - R exchange). These reservoir states populations monotonously increase at Tut < 1 and 
achieve maximal values at T* ~ T^ 1 ln(r/j/T). Then the state n population decreases as 

KWI 2 = ; 2_ L(r T l 7^ ( ex p(- K ) ~ cxp(-(r r + r,)t)) . (35) 

-K[n z + (T R + 7 S - T) z ) 
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• Regular Loschmidt echo regime, where stochastic-like behavior is not achievable, because the echo intensity is 
already negligibly small when the cycles start to mix 

i >r> 2^ ;r ^ 1>r - ^ 

In this regime we get Loschmidt echo phenomenon which occurs due to synchronisation of the backward (from 
the reservoir R into the system S) transitions. Furthermore when the state n of the reservoir is de-populated 
completely, its complex amplitude a n (t) acquires the phase shift ir. This leads to resonance behavior at the double 
reservoir R interlevel transition frequency. In turn these double transitions n — > s — > n' effectively redistribute 
the population between the reservoir R states. Besides the phase shift ir yields to a s (t) sign alternation between 
even and odd cycle numbers (see the Fig. 1). In odd cycle numbers the energy (or population) flows from the 
reservoir R into the system, whereas in even cycle numbers it flows from the system into the reservoir states. 

• Quasi-stochastic time evolution when 

when the system time evolution is described by incoherent dynamics with multicomponent Loschmidt echo. In 
this case the system exchanges its energy with many states of the reservoir. 



Cycle mixing regime 



r « — ; T R » 1 > T s . (38) 



In this case the cycles are strongly overlapped, the system amplitude oscillates in a shorter and shorter time 
scale, because the number of zeros of the a s (t)^ in the cycle k increases proportionally with k. Therefore to 
keep fully deterministic system description at t —> oo, one has to know the system evolution at shorter and 
shorter time intervals. This phenomenon leads to effective loss of deterministic dynamics at arbitrary small time 
coarse-graining, related physically to uncertainties of the system or of the reservoir, for example. 



We see from these conditions that the Loschmidt echo occurs depending mainly on the decay rate of the R states, 
whereas stochastic like behavior occurs depending on the product of T ■ T^, which characterises the energy (or level 
population) flow from the system into the environment Q through the reservoir R states as intermediate ones. We 
illustrate these findings in the Fig. 1 which shows how system level e s population (oc |a s (t)| 2 ) depends on the coupling 
to the reservoir R states. The S — Q coupling affects also the shape of the Loschmidt echo signal, suppressing upon 



Figure 1: System population time evolution for the model parameters: C 2 = 1, T s = 0.2. F shows how strong intramolecular 
vibrations are coupled to environment: (a) F — 0.01; (b) F = 0.05; (c) F = 0.15. Time is measured in the recurrence cycle 
period units 2tt/Q, where f2 is the interlevel spacing in the reservoir R. 

F s growing the most intensive echo component I^k ■ Our results demonstrate the main advantage of the partial cycle 
amplitude representation in comparison with the standard Fourier representation for the a s (t) evolution. The latter 
one depends on the complex manner of interplay of all three relevant coupling parameters S — Q, R — Q, and S — R, 
whereas in the former one these three interactions are factorised. 
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A. Averaged over cycle population of state 



The system S population averaged over a cycle k is determined naturally as 



271 



(|ai fc )| 2 )^ / dt'|a< fc >(t')| 2 . (39) 



o 



For non-overlapping cycles the upper integration limit in the (|39|) can be moved to infinity (since the integrand is 

exponentially small for t' > 2tt). Replacing ai k \t') in the (|39|) by its Laguerre polynomial explicit form (|27|) . and 
utilising known [4l| integrals with Laguerre polynomials we find 

<|a« | 2 ) fe * exp(-4,r)(l + -3^ + 1.3/2,-fc + 3/2,(l-C)/(l + C )) 3 ) 

fl vv A sy +1, 3/2, -fc + 3/2,1) ' v ' 

where F(a, b, c, z) is a hyper-geometric function [4l|, and £ = T s /Tr. When £ -C 1, the (|40|) becomes (^-independent, 
and the averaged over cycle system state population decays with the R — Q transitions rate constant 4wT. Let us 
stress again the message we got from the comparison of the expression (1241) for the initial cycle decay rate and the 
expression (|40[) for long-time evolution: the former is determined by the S — Q transitions, whereas the latter one by 
the R — Q transitions. 



B. Notes on the coarse grained spectrum 



As it was pointed out above, when k (cycle number) is growing, the number of the fine structure components of the 
Loschmidt echo signals grows proportionally to k. As a result of increasing k, the echo shape becomes more and more 
complex. To detect the Loschmidt echo signal accurately one has to use a detector with better and better resolution 
because the signal shape upon k growing includes more and more components. The same can be formulated that for 
a given detector resolution, there is always a certain cycle number, when the accuracy of the detector is not sufficient 
to resolve the Loschmidt echo signal. Therefore the finite accuracy of any detector device yields to the fact that 
one may not restore the energy spectrum by measuring the time evolution of the system amplitudes a s (t) . Although 
mathematically speaking to restore the spectrum it is enough to perform inverse Fourier transform of the amplitude, 
in practical terms it is not possible, since the amplitude cannot be measured with the required accuracy for high 
enough cycle number. With a finite detector resolution we are able to perform only truncated Fourier transformation 
up to a certain limiting recurrence cycle number threshold K m . We illustrate this message on a reduced accuracy of 
truncated Fourier transformation in the Fig. 2. Note that noise from the measurement device also restricts a spectrum 
part which can be reproduced from asymptotic behavior of a s (t). The Fig. 2c shows spectrum distortion induced 
by external white noise. Namely, we have calculated the R -reservoir absorption spectrum for the truncated at K m 




Figure 2: Adsorption spectrum for the reservoir R calculated by truncated Fourier transformation of the amplitudes for the 
model parameters F = T s — 0.1, C 2 = 1, and the limiting by measurement precision recurrence cycle number K m is: (a) 
K m — 1; (b) Km = 5; (c) K m = 20. In the Fig. 2c white noise produced by a measurement device (signal/noise ratio 10/1) is 
also added. 

Fourier transformation of the amplitudes. We see that coarse graining produced by spectral measurement device, 
reduces the reproducible size of the spectrum, and distorts strongly the line shape. In the Fig. 2a (K m is sufficiently 
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small) no traces from the recurrence cycles are visible, and we get Fermi Golden Rule exponential decay. For the 
larger values of the K m (Fig. 2b) there are clear indications for adsorption peaks corresponding to the Loschmidt echo 
signals in a s (t). However, adding external white noise (Fig. 2c) we see that the adsorption peaks are reduced even 
for the bigger (than in the Fig. 2b) values of K m . In its spirit the described mechanism is in the heart of quantum 
mechanical uncertainty produced by the interaction between quantum system and classical measurement device. It 
is considered as a route to quantum chaos and irreversibility phenomena [42| . 



VI. CONCLUSION 



To summarise, in this paper we investigated quantum dynamics of a single selected state of a small system. This 
selected state is coupled to other discrete dense states of the same system, and all the states are coupled also to 
continuous spectrum environment. This publication represents a substantial extension of the note (3l| . where only a 
single equidistant reservoir and constant-coupling Zwanzig model have been studied. Here we provide a more complete 
account and investigation of phenomena only briefly addressed in (3lj | , and generalise the Zwanzig model to include 
quasi-stationary states and coupling to the second continuous spectrum reservoir - environment Q. We show that 
there are possible four dynamic regimes of the evolution: 



• (i) - independent of the environment Q exponential decay suppressing backward R -S transitions 

• (ii) Loschmidt echo phenomenon occurring not only due to almost coherent oscillations governed by transitions 
from the system to the resonance reservoir R state, but also due to double transitions n — > s —* n' at the double 
reservoir R interlevel transition frequency 

• (iii) - incoherent dynamics with multicomponent Loschmidt echo, when the system is exchanged its energy with 
many states of the reservoir 

• (iv) - cycle mixing regime, when due to unavoidable coarse graining in any real system of time or energy 
measurements, or initial condition uncertainty, the system loses invariance with respect to time inversion. In 
such conditions dynamic evolution of the system cannot be determined uniquely from the spectrum, and in this 
sense long term system dynamics appear to be random. 



The quantum dynamics of the selected level demonstrates non-trivial fine structures of the recurrence cycles 
(Loschmidt echo) and cycle mixing leading eventually to irregular, chaotic-like long term evolution. Our results 
illustrate non-ergodic dynamics of such a system, i.e., system population (or its energy) is not equally distributed 
over all system states but in certain time intervals it is concentrated in a few levels. The generalised Zwanzig model 
investigated in this paper reflects the spirit of minimalist approaches, in that it is simple yet based on a physical 
principle. The results presented here are probably less notable in terms of technological applications of nano-systems, 
than for the progress they could generate in our understanding of their complicated vibrational spectra. 

Our consideration yields quite reasonable qualitative description of a variety of vibrational relaxation regimes 
and mode selectivity observed in experiments, and the model under investigation appears to be the simplest one 
demonstrating that relatively small variation of the coupling enables us to change qualitatively the dynamic behavior. 
One of the main difficulties in a quantitative comparison of the results of our simplified model with specific experimental 
measurements or elaborated numerical simulations is the availability of an accurate connection between experimental 
control parameters and entering theoretical model coefficients. In this sense our model should be treated as a working 
hypothesis. Although it provides a qualitative insight to the intra- molecular vibrational dynamics, we are left with 
many questions unanswered which must be perused in further work. Understanding all its limitations, we nevertheless 
hope that our theory captures the essential elements of intra-molecular vibrational relaxation in nano-systems. Note 



that modern femtosecond spectroscopy methods (see e.g., [17| - [23|, [43| - |44[) indeed demonstrate (in a qualitative 



agreement with our consideration) remarkably different types of behaviors (exponential decay and irregular damped 
oscillations) of relatively close in energy initially excited states. We believe that we are the first to explicitly address 
this issue. 
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